第 37 卷 第 5 期 
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小 波 系数 与 小 波 估计 系数 之 间 存 在 的 恒定 偏差 问 


题 ， 提 出 一 种 基于 改进 阀 值 函 数 的 图 像 去 噪 算 法 。 该 算法 结合 已 改进 阔 值 函数 的 优点 ， 通 过 设置 适当 的 调整 参数 动 


态 选 取 固 定 阅 值 ， 增 加 调节 因子 来 降低 原 小 波 系数 和 估计 ， 


小 波 系 数 之 间 存 在 的 恒定 偏差 ， 从 而 提高 重 构图 像 和 原 图 


像 的 逼近 程度 。 改 进 后 的 冰 值 函数 在 冰 值 处 满足 连续 性 ， 同 时 满足 函数 的 渐进 性 和 高 阶 可 导 性 。 仿 真 结果 表明 ， 采 
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Wavelet transform image denoising algorithm based on improved threshold function 


Zhang Huijuan, Zhang Damini, Yan Wei, Chen Zhongyun, Xin Ziyun 
(College of Big Data & Information Engineering, Guizhou University, Guiyang 550025, China) 


用 改进 后 的 阔 值 孙 数 进行 图 像 去 骂 , 视觉 效果 好 , 男 一 方面 对 均 方 误差 (MSE)、 峰 值 信 品 比 (PSNR) 和 信 骂 比 (SNR) 
等 去 品 性 角 0 PSNR 和 SNR 都 提高 了 ， MSE 有 所 降低 ， 去 噪 效果 得 到 了 优化 。 
; 峰值 信 骂 比 : 信 骂 比 


Abstract: Aiming at the discontinuity of the traditional hard threshold function at the threshold and the constant deviation 


between the original wavelet coefficient and the wavelet estimation coefficient in the soft threshold function, an image 


denoising algorithm based on the improved threshold function is proposed. The algorithm combines the advantages of the 


improved threshold function, dynamically selects the fixed threshold by setting appropriate adjustment parameters, and adds 


the adjustment factors to reduce the constant deviation between the original wavelet coefficient and the estimated wavelet 


coefficient, thereby improving the degree of approximation of the reconstructed image and the original image. The improved 
threshold function satisfies continuity at the threshold while satisfying the asymptotic and higher order conductibility of the 


function. Simulation results show that the improved threshold function is used for image denoising. There is a good visual 


effect. On the other hand, the denoising performance indicators such as mean Square error (MSE) , peak signal-to-noise ratio 
(PSNR) and signal-to-noise ratio (SNR) are performed. In comparison, both PSNR and SNR are improved, MSE is reduced. 


The denoising effect is optimized. 


Key words: Wavelet transform; threshold function; threshold image denoising; mean-square error; peak signal-to-noise 


Ratio; signal-to-noise Ratio 


0 引言 


了 一 定 的 改善 ; 文献 [9] 采 用 


妇 阔 值 函数 的 小 波 图 像 去 
因子 ， 但 3 


省 未 分 析 分 


提出 了 一 种 


去 噪 函数 的 影响 ， 文 献 


适应 的 阔 


噪 ， 根 
图 像 在 传输 的 过 程 中 不 可 避免 地 会 受 至 呈 度 的 噪 解 层 数 对 
污染 ， 造 成 区 问题 ， 影响 图 像 的 质量 吾 期 的 处 理 值 去 噪 函 数 ， 
等 。 因 此 ， 有 效 的 图 像 去 噪 方案 相当 重要 。 换 对 不 同 小 波 估计 系数 之 间 
频率 在 时 域 上 的 取样 长 度 可 以 进行 调节 ,具有 多 分 辨 的 特点 。 [11] 提 出 
1994 年 ,Donoho 也 nstone 在 小 波 变换 的 基 捉 出 了 小 设置 不 
波 阔 值 去 噪 算 法 , 该 方法 在 图 像 去 噪 中 的 应 日 渐 成 熟 上 4。 果 可 知 依然 存在 
型 的 的 硬 冰 值 函数 和 软 阔 值 函数 都 存在 相应 的 不 足 之 处 ， 闵 值 去 噪 函 
硬 阔 值 函数 在 闪 值 处 不 连续 ， 软 闵 值 函数 中 原 小 波 系数 和 小 真 结 果 可 知 ， 改 进 的 让 
波 估计 系数 之 | 间 存 在 恒定 的 偏差 间 题 。 针对 以 上 问题 ， 大 量 是 效果 不 是 特别 
学 者 对 闷 值 E EE 点 研究 ， 并 对 经 : 函数 进行 处 ， 本 文 综合 
了 改进 I 王 本 。 是 出 了 一 种 新 的 阔 值 函数 ， 弹 值 函数 值 函 数据 
中 存在 调节 医 虽然 可 以 动态 调整 阀 值 函 妆 去 噪声 并 比较 分 析 不 同 i 
效果 并 不 是 很 理想 ; td 种 改进 甬 数 ， 加 调节 因 时 
入 调整 参数 m 对 、 软 冰 值 函数 进行 优化 ， 效 果 得 到 分 解 层 数 ， 从 
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有 所 提高 ， 但 


收 稿 日 期 : 2018-10-09; 修 回 日 期 : 2018-12-03 


作者 简介 : 张 绘 娟 (1994-), 女 (土家 族 ) 
阳 人 ， 教 授 ， 博 士 ， 主 要 研究 方向 为 计 
硕士 研究 生 ， 主 要 研究 方向 为 计算 机 应 用 技术 、 优 化 计算 ; 陈 忠 云 (1989-)， 贵 州 遵义 人 ， 男 ， 
研究生， 主要 研究 方向 为 计算 机 应 用 技术 、 优 化 计算 . 


男 《通信 作者 )， 焉 州 贵 
男 ， 贵 州 贵阳 人 ， 
应 用 技术 、 优 化 计算 ; 辛 梓 其 〈1994-)， 贵 州 铜仁 人 ， 女 ， 硕 了 


算 机 应 用 技术 、 


献 所 提 方 案 存在 的 不 足 2 
函数 的 优点 和 已 改进 新 的 阔 


A 


设置 不 同 的 调节 因子 ， 


: 贵州 省 自然 科学 基金 项 目 〈 黔 科 合 基础 [2017] 1047 号 ) 

), 贵州 思 南 人 , 硕士 研究 生 , 主要 研究 方向 为 计算 机 应 用 技术 、 优 化 计算 、 信 号 与 信息 
智能 算法 、 信 号 与 信息 处 理 等 (1203813362@qq.com); 
硕士 研究 生 ， 主 要 研究 方向 为 计算 机 


层 数 来 选 l 
的 相关 调整 参数 和 最 佳 


再 选取 最 优 
色 最 佳 


` 处 理 ; 张 达 敏 《1967-)， 
闫 威 〈1993-)， 


201904.00006v1 


chinaXiv: 


录用 定稿 张 绘 娟 ， 
小 波 分 解 总 层 数 m。 


1 ”小 波 阐 值 去 噪 原理 及 其 计算 


1.1 小 波 阔 值 去 噪 原理 

小 波 阀 值 去 噪 算法 是 图 像 去 噪 中 常用 的 方法 ， 其 思想 
对 处 于 高 频 部 分 的 噪声 进行 闪 值 处 理 ， 再 对 处 理 后 的 小 波 系 
数 进 行 重 构 得 到 去 噪 后 的 图 像 。 阔 值 的 具体 估计 和 效 值 函数 
的 选取 是 影响 小 波 阔 值 函数 去 噪 效果 的 关键 因素 ， 阔 值 设 定 
太 小 ， 去 噪 后 的 图 像 仍然 存在 噪声 ， 反 之 太 大 则 会 将 重要 图 
像 特征 滤 掉 ， 造 成 重 构 图 像 模 糊 。 所 以 阐 值 函数 和 阐 值 判定 
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硬 冰 值 函数 定 义 如 下 : 
~ or Iwis| =4 
Wi 二 hw, < (4) 
软 阔 值 函数 定义 如 下 : 
sgnwi (wisl —A), wi =4 
on G) 
软 硬 闷 值 折 中 函数 定义 如 下 : 
人 sgn(w;y)(|wil —Q4), |w, =A, a=0.5 
Wi -| ho < (6) 
典型 硬 、 软 阔 值 函数 的 曲线 变化 如 图 2 所 示 , 从 图 可 知 ， 


准则 的 选取 尤为 关键 。 本 文 将 含 高 斯 白 噪 声 图 像 的 数学 模型 
描述 如 下 : 


f(j,D= gD +na(j,D) (1) 


若 采 用 硬 闹 值 函数 获取 的 小 波 估 计 系 数 连续 性 很 差 ， 小 波 估 
计 系 数 在 +6 处 是 不 连续 的 ， 因 此 最 后 重 构 得 到 的 图 像 会 产 
生 振 荡 ， 出 现 截 断 效 应 等 不 利 影响 ， 软 阔 值 函数 具有 较 好 的 


其 中 : 70 表示 含 高 斯 白 噪 声 的 图 像 ，g(7.0 表示 不 含 噪声 
的 图 像 ，n(j,D) 表示 高 斯 白 噪声 ， 服 从 正 态 分 布 N(0,57)，j 和 
/表示 图 像 的 像素 位 置 。 小 波 阔 值 去 噪 的 主要 步骤 (图 1) 如 下 : 

a) 对 A(j,D) 作 小 波 变换 ， 得 到 各 尺度 的 小 波 系数 wj, ; 

b) 设 定 各 分 解 层 靖 值 ， 对 小 波 系数 w 进行 闵 值 处 理 , 得 
到 小 波 估计 系数 如 v ; 
c) 重 构 小 波 系 数 , 再 利用 小 波 估计 系数 进行 重新 构造 


得 到 去 噪 处 理 后 的 图 像 广 。 


合 噪 | ,| 小波 闪 值 和 间 什 图 像 去 嗓 
图 像 | "| 分 解 | “| 函数 量化 处 理 | “| 重 构 图 像 | 


图 1 小 波 阔 值 去 噪 流程 图 
Fig. 1 Wavelet threshold denoising flow chart 

1.2 ”小波 阅 值 的 计算 

在 小 波 阔 值 去 噪 领域 中 ， 阔 值 的 选取 直接 影响 小 波 闵 值 
的 去 噪 效果 。 目 前 常用 的 阔 值 选取 法 有 : Stein 无 偏 似 然 估 计 
(rigrsure 准则 )、 启 发 式 阔 值 (heursure 准则 )、 最 大 最 小 六 
值 (minimaxi 准则 )、 固 定 闵 值 (sqtwolog 准则 )。 采 用 不 同 
的 阔 值 准则 会 产生 不 同 的 去 噪 效果 ， 固 定 阔 值 和 启发 式 阔 值 
在 确定 阔 值 时 ， 阔 值 选取 过 大 ， 导 致 图 像 的 细节 系数 容易 于 
失 ，Stein 无 偏 似 然 估 计 和 最 大 最 小 闷 值 在 确定 闵 值 时 ， 阔 值 
选取 过 小 ， 去 噪 效果 不 明显 。 目 前 常 使 用 的 是 Donoho 等 在 
高 斯 噪声 模型 基础 上 031， 通 过 独立 正 态 变量 决策 理论 而 得 出 ， 


连续 性 ， 但 是 小 波 系数 与 小 波 估计 系数 之 间 存在 一 定 的 偏差 
问题 ， 导 致 重 构图 像 与 原始 图 像 之 间 的 逼近 度 减 弱 ， 达 不 到 
理想 的 重 构 效果 ， 从 而 影响 最 后 的 去 噪 效果 。 


入 


个 
硬 阐 值 函数 


.8 一 ”~ 略 六 信函 数 
> 


pe Bb Ww 


图 2 人 硬 闵 值 函数 与 软 阐 值 函数 
Fig.2 Hard threshold function and soft threshold function 
2.2 ”改进 的 小 波 阐 值 函 数 
针对 传统 的 硬 软 阔 值 去 噪 函数 存在 的 缺陷 问题 , 文献 [17， 
18] 提 出 的 改进 阔 值 函数 分 别 如 下 所 示 : 


24 
a sgn(wi)x (|w,l 一 Iwis| =4 四 
0 Iwis| <4 
sgn(w;)x (|w,,l —A+A4/ (2k+1)), Iw 三 1 
Wi = Wi ww < (8) 
(Qk+D: A 


其 中 : 式 (7) 很 好 地 解决 了 含 噪 图 像 的 小 波 系数 和 估计 小 波 


其 表达 式 如 下 ， 
14.62 O) 


系数 存在 恒 等 差 的 问题 ， 但 是 阔 值 函数 中 不 含 调节 因子 ， 因 
此 不 能 动态 调整 闷 值 函数 ， 式 (8) 中 存在 调节 因子 kx， 可 以 
动态 调整 阅 值 函数 值 ，k 取 值 为 1。 并且 在 wi| <4 的 区 内 ， 


其 中 : 4 表示 闵 值 ，MxN 表示 图 像 的 大 小 ，5 表示 高 斯 噪声 
方差 ,实际 中 噪声 的 方差 是 无 法 预知 的 , 需要 事先 估计 得 到 ， 
通过 5=median(|wij)10.6745 来 获取 高 斯 噪声 方差 。 为 了 使 阀 值 
不 限定 在 某 一 个 固定 值 ， 本 文采 用 、 文 献 [14] 中 的 改进 阔 值 
选取 法 ， 其 数学 表达 式 如 下 : 


ev = (6V2In(M x N)) x eun) G) 


并 不 是 简单 地 将 闵 值 函数 置 为 0， 而 是 通过 一 个 非 线性 的 函 
数 实现 对 闵 值 的 逐渐 压缩 ， 既 可 以 保证 阀 值 函数 的 连续 性 又 
可 以 避免 阔 值 函数 的 直接 截断 而 引起 的 震 功效 应 ， 可 以 很 好 
地 解决 小 波 系数 与 估计 小 波 系 数 存在 的 恒定 偏差 问题 ， 具 有 
较 好 的 阔 值 去 品 效 果 。 

2.3 ”一 种 新 的 改进 小 波 阐 值 函数 

本 文 结合 传统 的 硬 、 软 阔 值 函数 的 优点 和 当前 已 改进 阔 


五 


下 文 涉及 到 的 固定 阔 值 都 由 式 (3) 来 更 新 ，ec 为 收缩 因 
子 ，n 为 小 波 总 分 解 层 数 。 
2 小波 阐 值 去 噪 函 数 的 改进 
2.1 常用 的 小 波 阐 值 函数 

阔 值 函数 的 选取 将 直接 决定 小 波 系数 的 处 理 策 略 ， 从 而 
影响 最 后 的 去 噪 效果 。 其 中 广泛 应 用 的 是 由 Donoho 等 提出 


值 函数 的 思想 ， 对 闵 值 函数 进行 改进 。 通 过 增加 调节 因子 来 
缩小 原 小 波 系数 与 估计 小 波 系 数 之 间 存 在 的 恒定 偏差 ， 因 此 
提出 一 种 新 的 改进 阔 值 去 噪 函 数 ， 其 表达 式 如 下 ; 


4 
uwis + (1—u) sgn(w;)(|w,,| pr )，|wj| >4 
WA w, 22 (9) 
sen 一 一 “一 一， V | <4 
sn( rep lw 


的 硬 、 软 阀 值 方法 ， 之 后 大 量 研究 学 者 在 此 基础 上 对 阔 值 函 
数 作出 改进 ， 例 如 软 硬 阔 值 折 中 函数 05 9， 传统 的 阔 值 函数 
表达 式 如 下 所 示 : 


其中 : =1-e20oH ，a，D，1， 为 可 调 参 数 ， 日 都 为 正 数 ， 
所 以 可 以 通过 选取 不 同 的 数值 来 提升 阅 值 函数 的 去 噪 性能。 
并 且 在 wal <4 的 区 间 内 ， 并 不 是 直接 将 阔 值 函数 置 为 0， 而 
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是 通过 一 个 非 线性 的 函数 实现 对 阅 值 的 逐渐 压缩 ， 这 样 可 以 


避免 阔 值 函数 的 直接 截断 而 引起 的 震荡 效应 ， 进 一 步 解决 阐 
值 函数 中 存在 的 恒定 偏差 问题 。 

从 数学 角度 去 考察 改进 的 阔 值 函数 : 

1) 函数 的 连续 性 分 析 


lim, lim [ uwis +(1—u) sgn(w)(|w;sl 一 
jh Wi 


)] 


Cr 


=,lim (1-0 sgnwi lw 一 Fa 


= ml sgn(w)(|w,,l Zap) 
4 


= lim [ w， 一 一 一 一 一 
il 及 Zo] 


=0 


lim w,,=0 
vj A+ J 


同 理 可 证 : J Wi,=0 6 


则 有 ,im 总 = ,lm 名 ,=0 ， 说 明 该 阔 值 函数 在 4 处 是 连 
续 的 ; 


,im, Wi =, lim, [uw + -0) Sen) ~ sap) 
= (lwal 一 TPR 
=, lim [一 (own| = 
=0 
a Wa=0, Sy w=0 


则 有 ,lm = lmeWs=0， 说 明 阔 值 函数 在 -4 处 也 是 


连续 的 。 

综 上 所 述 ， 改 进 的 阀 值 函 数 在 +4 处 都 是 是 连续 的 。 
2) 函数 的 渐 近 性 分 析 

当 w + 时， 


4 

Wn +(1—u) sgn(w;)(|w;,l a 
lim 一 一 = lim 
Wy Wi wi Wi 

4 
; Qo) (wl ~ ae) 

= lim [u+ ] 

Wj w, 

EE 
il 
当 Wt = = 时 » 

和 ， 村 4 
ey uw + (1—u) sgn(w;)(|w,l a ) 
lim = 
0 Wi Wj wi 

1 
. -1 (wil 一 PT) 
lim [zx 十 | 
Wm Wi) 
=1 
同时 ， 当 Wag > 时 : 
lim (WwW, —w,)) = 
Wo 
lim {[uw;, + (1—w) sgn(w;)(hw;il Sa E } 
Wi WW gn\W NW Wr i 


ut wis — Wir} 
a 
综 上 所 诉 ， 本 文 改进 的 阔 值 函 数 式 9) 是 以 和 =wii 为 
渐 近 线 的 。 
3) 改进 阐 值 函数 的 偏差 性 分 析 
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lim (Ww,,—w,,)= 
cn i 1) 


lim low + 4-0) sgnw)lwnsl ~ za ) Wi] 
= lim [wj,—w,] 

Wo 
=0 

lim (Wj—w;)= 
Wj 
im, low +(1—u) sgn(w)(|ws| — DR ) 一 whi] 
= lim [wj,—w,] 

Wj- 


=0 


所 以 可 得 lim wi-wi)=,lim (wn-wi)=0，, 当 wj, 逐渐 增 大 


时 ， Wi 和 wi 之 间 的 偏差 会 逐渐 减 小 ， 这 样 可 以 克服 软 闵 值 
函数 中 存在 的 固定 偏差 问题 。 
4) 函数 的 高 阶 可 导 性 
从 改进 的 阔 值 函数 的 表达 式 可 知 ， 当 wl>4 时 ， 改 进 的 
闷 值 函数 满足 高 阶 可 导 ， 有 利于 进行 各 种 数学 变换 及 处 理 。 
5) 阔 值 可 调 因 子 a，b，m，t 的 影响 分 析 
当 Mwl>4 时， 当 a=0 且 b=0 时 ,改进 阐 值 函数 为 软 闪 值 
函数 ， 当 2 一 时 ， 改 进 阔 值 函数 为 硬 阔 值 函数 。 因 此 ， 改 
进 的 阔 值 函数 可 以 在 软 、 硬 阔 值 函数 之 间 进 行 调整 ， 可 以 根 
据 实 际 应 用 来 增 大 或 者 减 小 a, bp»， 的 值 ， 进 行 合理 调节 ， 从 
而 达到 最 佳 去 噪 效果 。 在 wi <4 的 区 间 内 ， 并 不 是 直接 将 阐 
值 函数 置 为 0， 而 是 采用 一 个 非 线 性 的 函数 对 闪 值 进行 逐渐 
压缩 ， 同 时 通过 选取 不 同 的 上 值 可 以 调整 阔 值 的 取 值 。 为 了 
获取 更 优 的 调整 参数 ， 根 据 经 验 取 值 可 知 ，a,，b，m 都 为 正 
数 ，te[0 3] 为 正 整数 ， 并 且 通 过 反复 实验 比较 得 到 ， 当 可 调 
参数 a 取 值 为 1, b 取 值 为 1，m 取 值 为 2, t 取 值 为 1 时, 去 
噪 效果 最 佳 。 
综 上 ， 通 过 从 改进 闵 值 去 噪 函 数 的 连续 性 、 渐 近 性 、 仿 
差 性 以 及 调节 参数 进行 分 析 ， 从 给 出 的 证 明 可 知 ， 本 文 改进 
的 阔 值 函数 进一步 殉 服 了 硬 阔 值 函数 存在 的 不 连续 性 和 软 阔 
值 函数 中 小 波 系数 存在 的 恒定 偏差 问题 。 
基于 改进 小 波 阔 值 函数 的 图 像 去 噪 算法 步骤 如 下 : 
a) 选择 合适 的 小 波 基 函 数 db5， 然 后 对 含 噪声 的 图 像 
f(j) 进行 多 层次 的 正 交 小 波 变换 ， 获 取 一 组 小 波 分 解 系数 


Ww 


il; 


b) 根 据 5=mediandwj)10.6745 对 噪声 方差 5 进行 估计 ， 在 
此 基础 上 ,使 用 式 (3) 的 定义 来 计算 固定 阔 值 4 ， 以 4 为 门 
限 值 ， 通 过 式 〈9 ) 对 小 波 分 解 系数 进行 阔 值 处 理 ， 然 后 获取 
小 波 估计 系数 六; 
c) 根 据 阔 值 处 理 后 获取 估计 小 波 估计 系数 名 和 处 于 低 
频 部 分 的 小 波 系数 , 然后 利用 小 波 基 函 数 来 进行 小 波 重 构 后 ， 


得 到 进行 阔 值 去 噪 后 的 图 像 记 .本文 改进 算法 的 流程 图 如 图 


3 所 示 。 
当 |w| 24 和 |wi| <4 时 ， 令 表达 式 分 别 如 下 所 示 : 


Is 


) 


和 
Win =uwj, +(l—u) sgn(w;, wl 一 Capa 


2f+1 
Wi 


(e: 一 1/2)42 


3 ”实验 结果 及 其 性 能 分 析 

为 了 验证 本 文 提 出 的 改进 阔 值 去 噪 函 数 的 有 效 性 ， 采 用 
小 波 基 函数 db5 进行 不 同 层 数 的 分 解 与 重 构 。 所 选用 的 测试 
图 像 大 小 都 是 512x512 。 对 含有 不 同方 差 值 噪声 的 Lena 图 像 


人 
wi2 = Sgn(w;i,) 
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和 peppers 图 像 ， 分 别 使 用 典型 的 硬 、 软 阔 值 函数 、 软 硬 折 
中 准 值 函数 、 文 献 [17，18] 和 本 文 改进 的 闵 值 函数 进行 实验 
测试 。 选 用 不 同 的 高 斯 噪声 方差 值 : 5=0.01, 5=0.03 , 5=0.05， 
对 图 像 进行 加 噪 ,采用 小 波 基 函 数 db5 进行 小 波 分 解 与 重 构 ， 
同时 选用 不 同 的 分 解 层次 数 n。 当 5=001，n=3 时 ，Lena 去 
图 像 和 Peppers 去 噪 图 像 如 图 4、5 所 示 ， 从 视觉 角度 上 可 
以 看 出 ， 本 文 改 进 阔 值 函数 下 的 去 噪 图 像 比 另 外 五 种 去 噪 函 
数 下 的 去 噪 图 像 更 清晰 ， 去 噪 效果 明显 。 

为 了 进一步 验证 本 文 所 提出 的 改进 阔 值 函数 的 去 噪 效果 ， 
从 客观 数据 上 比较 分 析 所 提出 的 阔 值 函数 的 性 能 。 使 用 均 方 
误差 (mean squared error，MSE )、 峰 值 信 噪 比 (peak 
signal-to-noise ratio，PSNR ) 609 和信 噪 比 CSNR) P20 对 去 噪 
后 的 图 像 进行 数值 计算 , MSE 的 数值 越 小 说 明 去 噪 图 像 质量 
越 好 ;峰值 信 品 比 的 数值 越 大 ， 说 明 去 品 效 果 越 好 。 均 方 误 
差 的 表达 式 如 下 所 示 : 


源 


1 刀 嫌 ,ra 
MSE= x [OD 7) ] (10) 
峰值 信 噪 比 和 信 噪 比 的 表达 式 分 别 为 
255? 
PSNR =10x le(sE) (11) 
SNR=10xlg{[2, 2 7 C.D)°]/ M4 x N)- MSE} (12) 


其 中 : MxN 为 测试 图 像 大 小 ，7(jD) 和 ,分别 表示 为 无 品 


图 像 和 去 噪 图 像 。 为 了 更 全 面 地 比较 阔 值 函数 的 性 能 ， 进 
步 比 较 分 析 不 同 去 品 方 法 下 的 均 方 误差 (MSE)、 峰 值 信 噶 


像 ”(c) 硬 闵 值 去 噪 ”(d) 软 阐 值 去 噪 


比 (PSNR) 和 信 噪 比 (SNR )。 选 月 
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目 了 不 同 的 噪声 方差 5 和 


不 同 的 分 解 层次 数 n， 通 过 实验 仿真 得 到 的 数据 如 下 图 表 所 


不 。 


开始 


大 


含 噪 图 像 f(j,1 ) 进行 多 
层次 的 正 交 小 波 变换 


公式 〈3) 计算 固定 阅 值 4 


通过 公式 〈9) 对 小 波 分 解 系数 
行 阅 值 处 理 


通过 对 阅 值 处 理 后 得 到 的 小 波 估计 系数 和 
处 于 低频 的 小 波 系 数 ， 利 用 小 波 基 函数 


db5 进 行 小 波 重 构 


Fig.3 


(e) 软 硬 折 中 去 品 ”(D 文 献 [17] 去 吕 


(a) 原 始 图 像 ” (b) 含 噪声 图 像 ”(c) 硬 闵 值 去 品 ”(d) 软 闵 值 去 噶 


图 $ Peppers 图 像 采 


采用 不 同 的 阔 值 去 噪 方法 对 加 噪 的 Lena 图 像 进行 阔 值 


图 4 Lena 图 像 采用 不 同方 法 去 噪 效 果 


Fig.4 Lena image uses different methods to denoise the effect 


(e) 软 硬 折 中 去 噪 (文献 [17] 去 噪 


用 不 同方 法 去 噪 效果 


输出 去 噪 后 的 图 像 


图 3 改进 算法 流程 图 


Improved algorithm flow chart 


Fig.5 Peppers image uses different methods to denoise the effect 


(g) 文 献 [18] 去 噪 


(hb) 本 文 去 噪 


(hb) 本 文 去 噪 


噪声 方差 值 增加 到 0.11 时 ，6 种 阔 值 去 噪 函 数 下 的 PSNR 都 


去 噪 ， 评 定 指标 数值 如 表 1 所 示 。 从 表 中 数据 的 变化 趋势 可 
知 ， 与 传统 的 硬 软 阔 值 去 噪 函 数 相 比 ， 当 前 已 改进 的 冰 值 去 
噪 函 数 下 加 噪 图 像 的 峰值 信 噪 比 PSNR 增加 了 约 8 dB, 均 方 
误差 MSE 值 大 大 降低 了 ， 信 噪 比 SNR 增加 了 约 7.5 dB 。 本 
文 改 进 阀 值 函数 下 的 峰值 信 噪 比值 最 大 值 为 30.734 9， 与 传 
统 的 软 阔 值 函数 下 的 PSNR 相 比 ， 增 加 了 约 11 dB， 由 此 可 
知 ， 采 用 本 文 阔 值 函 数 下 的 去 噪 效果 是 最 佳 的 。 

同 理 ， 采 用 不 同 的 阐 值 去 吕方 法 对 加 品 peppers 图 像 进 
行 闵 值 去 品 ， 去 噪 图 像 的 评定 指标 数值 如 表 2 所 示 。 

针对 Peppers 加 噪 图 像 ， 取 同一 小 波 总 分 解 层 数 n=3, 设 
置 不 同 的 噪声 方差 ， 峰 值 信 噪 比 PSNR 随 着 噪声 方差 值 的 变 
化 而 变化 ， 仿 真 结 果 如 图 6 所 示 。 从 曲线 的 变化 趋势 可 知 ， 
随 着 噪声 方差 的 增 大 ， 峰 值 信 噪 比值 PSNR 在 逐渐 降低 。 当 


一 定 程度 地 减 小 了 ， 与 其 他 阔 值 去 噪 函 数 下 的 PSNR 相 比 ， 
本 文 闪 值 函数 下 的 PSNR 仍然 是 最 大 的 ,去 噪 效果 比较 理想 。 


表 1 舍 品 Lena 


像 去 噪 效果 评定 指标 (n=3, 5=0.01 ) 


Table 1 Noise-free Lena image denoising effect evaluation index (n=3, 


5=001) 

去 噪 方法 MSE PSNR/dB SNR/dB 
含 噪声 图 像 640.5461 20.0653 14.4081 
硬 疮 值 去 噪 138.6119 26.7128 21.0557 
软 阔 值 去 噪 182.4737 25.5188 19.8616 

软 硬 折 中 去 品 157.1739 26.1670 20.5099 
文献 [17] 去 噪 105.4922 27.8986 22.2414 
文献 [18] 去 噪 114.5066 27.5425 21.8853 
本 文 阔 值 去 噪 54.9023 30.7349 24.6808 
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表 2 含 噪 Peppers 


图 


像 去 噪 效果 评 


张 绘 娟 


定 指标 (n=3, 6=0.01 ) 


Table 2 Noisy peppers Image Denoising Effect Evaluation Index (n=3, 


56=0.01 ) 
去 噪 方法 MSE PSNR/dB SNR/dB 
含 噪声 图 像 626.6721 20.1604 14.3506 
硬 阔 值 去 噪 121.5823 27.2821 21.4724 
软 阔 值 去 噪 151.0220 26.3404 20.5306 
软 硬 折 中 去 噪 112.8523 27.6057 21.7959 
文献 [17] 去 噪 103.0980 27.9983 22.1885 
文献 [18] 去 噪 107.0065 27.8367 22.0269 
本 文 闷 值 去 噪 45.2533 31.5743 24.0783 


PSNR/dB 


一 


图 6 含 噪 peppers 


方差 


20 1 1 1 1 1 1 是 1 D> 
0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 01 0.11 


像 的 PSNR 随 噪声 方差 的 变化 


线 (n=3) 


Fig.6 PSNR ofnoise-containing peppers image with noise variance 


curve (n=3) 


同 理 ， 针 对 同一 peppers 加 噪 图 像 ， 当 小 波 总 分 解 层 数 
n=3 时 ， 取 不 同 的 噪声 方差 值 ， 信 噪 比 SNR 随 着 噪声 方差 值 
的 增 大 而 减 小 ， 仿 真 结果 如 图 7 所 示 。 从 曲线 的 变换 趋势 可 
知 ， 当 噪声 方差 增 大 时 ，6 种 阔 值 去 噪 函数 下 的 SNR 都 相应 


地 降低 了 ， 当 方差 为 0.11 时 ， 


比 另 外 5 种 方案 下 的 SNR 都 大 。 


图 


知 ， 


波 总 
图 像 中 


本 文 改 进 冰 值 函数 下 的 SNR 


一 合 一 本文 闵 值 去 噪 | .| 


7 今 


针对 同一 Peppers 加 噪 图 
同 的 小 波 分 解 层次 数 n， 
的 变化 而 变化 ， 


图 


含 噪 peppers 


访 真 结 呈 


像 的 
Fig.7 SNR ofnoise-conta 


curve (n=3) 


果 如 图 8 所 示 。 


小 波 分 解 总 
分 解 层 数 的 


Ml 


数 不 同 RT 
大 ， 峰 值 


的 有 用 信息 


本 文 改进 
的 PSNR 都 大 ， 去 


将 会 丢失 。 


理想 。 


涯 值 信 噪 比 PSNR 


到 的 PSNR. 


15 | | | | | | | | | 
0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 
方差 


SNR 随 噪声 方差 的 变化 曲 


线 (n=3) 


ining Peppers image with noise variance 


像 ， 取 噪声 方差 5-001 ， 设 置 不 
随 着 总 分 解 层 数 n 


从 


线 的 变化 趋势 可 
值 也 不 同 ， 随 着 小 
言 品 比 PSNR 呈 下 降 趋 势 ， 加 噪 
当 取 相同 的 小 波 分 解 层 数 n 时 ， 


的 阔 值 函数 下 的 PSNR 比 另 外 五 种 阔 值 去 噪 函数 下 
噪 效果 也 比较 
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PSNR/dB 


一 人 一 硬 阐 值 


图 8 


含 品 peppers 


6 


( 6=0. 


8 
小 波 总 分 


9 
小 解 层 数 


01 ) 


10 11 


加 像 的 PSNR 随 小 波 总 分 解 层 数 的 变化 趋势 


Fig.8 PSNR of the noisy Peppers image with the total number of 


wavelet decomposition layers ( 6=0.01 ) 


里 ， 针 对 同一 peppers 加 噪 图 


同 


像 ， 取 噪声 方差 5=001， 


设置 不 同 的 小 波 分 解 层 次 数 na， 得 到 的 信 噪 比 SNR 也 不 同 。 


仿真 结果 如 图 9 所 示 ， 


解 总 层 数 的 增 大 ，SNR 数 


从 


值 呈 下 降 趋 势 ， 当 


层 数 n 
函数 下 的 SNR 都 大 ， 


时 , 本 文 改 进 阔 值 函 


线 的 变化 趋势 可 知 ， 随 着 小 波 分 
取 同 一 小 波 分 解 


数 下 的 SNR 比 另外 5 种 阔 值 去 噪 
说 明 去 噪 效 果 比 较 理 想 。 


SNR/dB 


一 < 一 硬 阅 值 


时 9 含 噪 peppers 图 


wavelet decomposition layers ( 6=0.01 ) 
的 噪声 方差 和 不 同 的 小 波 总 分 解 层 
实验 数据 比较 分 析 


综 上 所 述 ， 设 置 不 同 
次 数 会 带 来 不 同 的 去 噪 效果 。 从 以 上 仿 
可 知 , 峰值 信 噪 比 PSNR 和 信 噪 比 SNR 随 着 噪声 方差 值 的 增 
有 分解 层 数 n 去 噪 效 果 不 一 样 。 
知 ， 本 文 所 提 的 阀 值 去 噪 函 数 在 一 
言 品 比 ， 降 低 了 均 方 误差 值 ， 去 


大 而 减 小 ， 选 取 不 同日 


5 种 阐 值 去 噪 函数 比较 可 


出 了 一 种 改进 上 
同 的 噪声 方差 和 不 同 的 小 波 分 解 层 数 ， 
去 噪 函 数 、 软 硬 折 中 浆 值 去 噪 函 


8 9 

小 波 总 分 解 层 数 

像 的 SNR 随 小 波 总 分 解 
( 6=0.01 ) 
Fig.9 SNR ofthe noisy peppers image with the total number of 


定 程度 上 提高 了 图 像 的 峰值 


噪 效果 有 明显 的 提升 。 
4 ”结束 语 
本 文 提 


层 数 的 变化 趋势 


与 男 外 


的 闹 值 图 像 去 噪 函数 ， 通 过 选取 不 


对 典 


出 的 阔 值 去 品 函 数 下 的 三 


个 评定 


闷 值 取 值 采 


用 了 新 的 闹 值 
可 以 选择 不 同 的 闵 值 ， 增 加 了 阐 


区 值 规则 ， 


型 的 软 、 硬 阔 值 
数 、 文 献 [17，18] 和 本 文 提 
指标 进行 比较 分 析 。 本 文 的 
对 此 ， 不 同 的 分 解 层 数 


值 选取 的 灵活 性 。 


文 改 进 的 阔 值 去 噪 函 数 含有 


数 ， 获 取 最 佳 的 小 波 估计 系数 。 


个 请 


视觉 效果 还 是 从 三 


调节 


闪 值 


FE 定 指标 


进行 


同时 ， 本 


刀子 ， 通 过 改变 调节 参 


仿真 实验 结果 


表明 ， 不 论 从 
比较 分 析 ， 在 去 除 噪声 


录用 定稿 张 绘 娟 ， 等 : 
的 同时 也 有 效 地 保留 了 有 用 信号 ， 与 其 他 去 噪 方案 相 比 ， 基 
于 本 文 改 进 阔 值 去 噪 函数 下 的 去 噪 效果 更 加 理想 。 
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